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Abstract 

We present the first femtosecond studies of electron-phonon (e-ph) thermalization in heavy 
fermion compounds. The e-ph thermalization time r ep increases below the Kondo temperature by 
more than two orders of magnitude as T = K is approached. Analysis using the two-temperature 
model and numerical simulations based on Boltzmann's equations suggest that this anomalous 
slowing down of the e-ph thermalization derives from the large electronic specific heat and the 
suppresion of scattering between heavy electrons and phonons. 
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Recent experiments have demonstrated that femtosecond time-resolved optical spec- 
troscopy is a sensitive tool to probe the low energy electronic structure of strongly cor- 
related electron systems [1-4], complementing conventional time-averaged frequency-domain 
methods. In these experiments, a femtosecond laser pulse excites a non-thermal electron dis- 
tribution. This non-thermal distribution rapidly thermalizes through electron-electron (e-e) 
interactions resulting in a change in the occupied density of states (DOS) in proximity to the 
Fermi energy (E/). Therefore, by measuring photoinduced reflectivity or transmissivity dy- 
namics as a function of temperature (T), it is possible to sensitively probe the nature of the 
electronic ground state. For example, femtosecond measurements of the carrier relaxation 
dynamics of high-T c superconductors and charge density wave compounds have provided 
new insights into the low energy electronic structure of these materials [1-3]. What is par- 
ticularly important is that even though the probe photon wavelength in these experiments 
ranges from the far-infrared Q, the near-I R Q, up to several evfl , the relaxation dynamics 
on identical samples is the same B, S npp„ rt ~ g the b a S1 c ld eaQ that dynamic photoinduced 
reflectivity measurements, in many instances, probe relaxation and recombination processes 
of quasiparticles in the vicinity of E/. 

In this Letter, we present the first studies of carrier relaxation dynamics in heavy fermion 
(HF) compounds using femtosecond time-resolved optical spectroscopy, aiming to elucidate 
the effect of localized /-electrons |5] on the quasiparticle relaxation dynamics. We have 
measured the time-resolved photoinduced reflectivity AR/R dynamics as a function of tem- 
perature on the series of HF compounds YbXCu4 (X=Ag, Cd, In) Jf| in comparison to their 
non- magnetic counterparts LuXCu4. Our results reveal that the carrier relaxation dynamics 
are extremely sensitive to the low energy DOS near Ep. In particular, in HF compounds 
the relaxation time r r shows a hundred-fold increase between the Kondo temperature (Tk) 
and 10 K, while in the non-magnetic analogues r r is nearly constant, similar to conventional 
metals like Ag, Au, and Cup]. Our analysis shows that the relaxation dynamics can be 
attributed to e-ph thermalization, and that the anomalous slowing down of the e-ph ther- 
malization stems from the large electronic specific heat in HF compounds and suppression 
of e-ph scattering within the peak in the enhanced density of states near Ey. 

In the following, we focus on YbAgCu 4 (a prototypical HF system with T K ~ 100 K 
and low temperature Sommerfeld coefficient 7 ~ 210 mJ/mol K 2 )[(| compared to its non- 
magnetic counterpart LuAgCu4 (7 ~ 10 mJ/mol K 2 ). The experiments were performed 
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on freshly polished flux-grown single crystalsp. We used a standard pump-probe set-up 
with a mode-locked Ti: Sapphire laser producing 20 fs pulses centered at 800 nm (photon 
energy hw p h — 1.5 eV) with an 80 MHz repetition rate. The photoinduced (PI) changes 
in reflectivity AR/R were measured using a photodiode and lock-in detection. The pum 
fluence was kept below 0.1 //J/cm 2 to minimize the overall heating of the illuminated spot 
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while the pump/probe intensity ratio was ~ 30. Steady-state heating effects were accounted 
for as described in [11], yielding an uncertainty in temperature of ±3 K (in all the data the 
temperature increase of the illuminated spot has been accounted for). 

Figure 1 presents the PI reflectivity traces on the two compounds at several temperatures 
between m 10 K and 300 K. The relaxation dynamics of the non-HF compound LuAgCu4 
display a very weak temperature dependence, with AR/R recovering on a sub-picosecond 
timescale at all T. The dynamics are similar to regular metals such as Au and Ag [3], where 
the recovery is predominantly due to e-ph thermalization. In contrast, Fig. 1(b) shows 
that for YbAgCu, the quasiparticle dynamics are strongly T-dependent. Specifically, above 
~ 140 K, the recovery time r r (determined by a exp(— t/r r ) fit to the data) is virtually 
T-independent but increases by more than two orders of magnitude as T — > K. We have 
measured similar dynamics on YbCdCu4 (TV ~ 100 K, 7 ~ 200 mJ/mol K 2 ). Furthermore, 
a similar divergence of r r occurs for CeCoIns below ~ 60 K, implying that the observed 
increase in the relaxation time starting at ~ TV and its subsequent divergence as T — > K 
is a generic feature of HF compounds and derives from their low energy electronic structure. 

The rise-time dynamics are also different in the two compounds. For LuAgCu4, the rise- 
time is ~ 100 fs at all temperatures. This is again similar to what has been measured 
on conventional metals and reflects the time it takes for the initially created high energy 
quasiparticles to thermalize towards Ef. Above ~ 25 K, YbAgCu4 displays a similar (fast) 
rise-time. Below this temperature the rise-time increases and, as the semi-log plot in Fig. 
1(b) reveals, becomes two-exponential at the lowest temperatures. Similar behavior also 
occurs for CeCoIn 5 , but is absent in YbCdCu4 indicating a strong dependence on the de- 
tails of the low energy electronic structure in HF compounds. While noting the presence 
of these anomalous rise-time dyna mics, further systematic studies are needed to obtain a 
more complete understanding. 12] In the following, we focus on the anomalous tempera- 
ture dependence of recovery dynamics below T^ which seem to be a general feature of HF 
compounds. 
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In conventional metals, the initial photoinduced change in the reflectivity arises from 
changes in occupation near Ej after e-e thermalization. Subsequently, the PI reflectivity 
recovery dynamics proceed on a picosecond timescale governed by e-ph thermalization |7[. 
The two temperature model (TTM) serves as a useful starting point in describing e-ph 
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thermalization in metals. [3 13, HQJ. The TTM describes the time evolution of the electron 
(T e ) and lattice (TJ) temperatures by two coupled differential equations^, In the low 
photoexcitation energy density limit, as in our case, when T e — T\ <C TJ over the entire 
temperature range, the set of two coupled differential equations can be linearized resulting 
in the following expression for the e-ph thermalization time 
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g{C: l + C^) . (1) 



Here C e and C\ are the electronic and lattice specific heats, respectively, and g(T{) is the 
e-ph coupling function j?l Oj- I n the case of simple metals, when the electron bandwidth is 
much larger than the Debye energy hw^, = /c^B^, and using the Debye model for the e-ph 
interaction, g(T) has particularly simple form. It is given by g(T) = dG(T)/dT, where 

/ rp \ 5 r e D /T 4 

G(T) = 4 9oo i^—j dx— T x(x,T). (2) 

Here g^ is termed the e-ph coupling constant, while x( x , T) is included to account for the 
variation in the electronic DOS, D e (e), and the normalized e-ph scattering strength F(e, e'), 
over the energy range Ef ± FujJd. It can be shown using Fermi's golden rule that 
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X(x, T) = - J de ew e D y ^ {/ (e) - f (e')} , (3) 
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where e' = e + ^, and ^ = xT and fo is the Fermi-Dirac distribution. In metals like Au or 
Ag, D d (e) and F are approximately constant in this energy range, i.e. D e (e) = D and 
F = 1, giving x = 1- fl'oo is typically 10 15 - 10 16 W/mol K (e.g. for Cu g^ = 6.2 x 10 15 
W/mol K corresponding to #(300fr) = 1 x 10 17 W/m 3 K Q)- 

At high temperatures (T > ©d), r ep (T) given by the TTM has been found to describe 
the temperature as well as photoexcitation intensity dependence of measured r r (T) j?l Q ■ 
Moreover, since the absolute value of r ep is determined by a single parameter g^, the 
technique has been successfully used to determine the dimensionless e-ph coupling constants 
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14j . However, at low temperatures (T < ©d/5) the TTM prediction 



of T ep oc T~ 3 has never been observed in metals 7] - instead, at low temperatures r r saturates 
at a constant value. The discrepancy between the experimental results and the TTM was 
found to be due to the fact that the TTM neglects e-e thermalization processes (by implicitly 
assuming that a Fermi-Dirac distribution is created instantly following photoexcitation). 
From simulations using coupled Boltzmann equations, Groeneveld et al. 7] showed that this 
discrepancy is due to the fact that at low temperatures the e-e and e-ph thermalization 
times are comparable. Since r ep oc T above ~ ©d/6 3,13, while r ee a T~ 2 - see Eq.(16) of 
Ref.Ql, the TTM is expected to fail at low temperatures where r ee <L r ep . 

In Fig. 2 we plot the T-dependence of r r on LuAgCu4 (solid circles), together with the 
TTM prediction for r ep (dashed line) given by Eqs.(P2D wi *h 6 D = 280 Kp, measured 
C e (T) and Cj(T) - see inset to Fig. 2, and = 2.6 x 10 15 W/mol K. Similar to Au 
or Ag[7], we find good agreement at T > 200 K, while below 40 K instead of showing a 
r ep oc T~ 3 divergence, r r saturates. 

In order to explain the discrepancy, we have carried out numerical simulations using cou- 
pled Boltzmann equations[3,Q|. Here, for example, the net phonon absorption by electrons 
with energy e is represented by [^f]^* = / dujK ep S(e, uj)D p {uj) D e (e + uj) , where D e (e) 
and D p {uj) are the electron and phonon DOS, and S(e,ui) = / e +w(l — fe) — ^w(/e — fe+u), 
with / and b being the electron and phonon distribution functions, respectively. K ep in 
the above equation and K ee in e-e scattering represent the square of the scattering ma- 
trix element, combined with all other numerical factors 3, When performing numerical 
simulations, a thermal phonon distribution (b t=0 = 6 (^1)) an d a non-thermal electron dis- 
tribution (f t =o = fo{T e ) ± 5f) was taken as the initial condition just after the laser pulse 
while r ep is found by fitting the total electron energy versus time curve to an exponential 
decay function. The initial perturbation 5f is around 10 -5 ~ 10~ 3 for the energy range 
between 0.10 ~ 0.15 eV above and below Ep, which is small enough that the increase in the 
temperature after e-ph thermalization is less than 1 K over the whole temperature range - 
consistent with the small excitation intensity used in the experiment fl^j] . The phonon and 
electron DOS used in the simulation were chosen such that they fit the specific heat data 
(i.e. for the phonon DOS, we use the Debye model D p (u) ~ uj 2 with hjj£, =24 meV, while 
D e (Ep) w 2.1 eV -1 f.u. _1 spin _1 ). The result of the simulation using the absolute value of 
K ep = 0.93 ps _1 eV and K ee /K ep = 700 is plotted by open circles jlTJ in Fig. 2. As expected, 
the simulation gives the same result as the TTM at high-T, while at low-T r ep saturates in 
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agreement with the experimental r r . 

Figure 3 shows r r (T) obtained on YbAgCu4. At high temperatures (T > Tk) the value 
of r r is similar to LuAgCu4. At low temperatures, however, r r increases by more than 
2 orders of magnitude. Since heavy fermions are characterized by a peak in the DOS at 
E'/P, the appropriate D el (e) should be used when modeling r ep (T). In our calculation 
we used D e (e) = -D pea fcexp[— (e/A) 2 ] + D , where D pea k = 70 eV _1 f.u. _1 spin _1 , A = 13 
meV and D = 2.1 eV _1 f.u. _1 spin _1 (identical to the value for LuAgCu^. It reproduces the 
experimental T-dependence of C e , as shown in the inset to Fig. 3. For simplicity we choose Ef 
at the center of the peak, so that the chemical potential is constant. Since D e (Ef) is almost 
two orders of magnitude larger than in LuAgC^ we expect that the e-e thermalization is 
much faster in YbAgCu4, and that the TTM would be valid at the lowest temperatures. 

The calculated r ep (T) using Eq.(0) is plotted in Fig. 3 by the dashed line. Here the 
approximate C e (T) and C\ (T) were used, g^ was taken to be the same as for LuAgCu 4 , 
while x ( X ,T) was evaluated explicitly for the above D e (e) and F — 1. Since r~ p a D e , and 
D e (Ef) ^> Dq the result is not surprising, implying that the simple TTM cannot account 
for the observed dramatic increase in r r at low temperatures. We should note that neither 
the value of the e-ph coupling constant nor D , which determine the absolute value of 
r ep , are necessarily the same in YbXCu4 and LuXCu4. However, even if the e-ph coupling 
is 10 times smaller in YbAgCu 4 compared to LuAgC^ (which would give 10 times larger 
value of r ep - as plotted by dashed line in Figure 3), the observed T-dependence of r r still 
cannot be accounted for. 

In order to account for the observed r r (T) we have to consider the nature of the electronic 
states within the peak in the DOS. In heavy fermions the peak in DJe) originates from 
hybridization of the localized /-levels with the conduction band electrons j5|. We hypothesize 
that the e-ph scattering within the DOS peak is suppressed, since the band dispersion near 
Ef is much weaker than in regular metals. Therefore it is quite possible that the Fermi 
velocity v p is smaller than the sound velocity v s , in which case momentum and energy 
conservation prohibit e-ph scattering when both initial and final electron states lie within 
the energy range where Vp < y a . m Assuming a parabolic band with Ep ~ T# ~ 100 K, and 
0.85 carriers per formula unit 18(, one obtains Vp ~ 4 km/sec, while the longitudinal sound 
velocity for YbIn 1 _ x Ag a .Cu4 (x < 0.3) is ~ 4.4 km/sec along [111] direction[19] (similar 
v s is expected for YbAgCu4). Even though a parabolic dispersion relation is just a rough 
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approximation, and a direct measurement such as de Haas-van Alphen effect is required to 
obtain vp, our simple estimate supports this idea. 

Using this hypothesis, good agreement with the data can be obtained. r ep (T) obtained 
by numerical simulations based on Boltzmann equations with K ep set to for processes 
where the initial and final electronic state are in the range of —24 < e < 24 meV (i.e. within 
the DOS peak), and K ep = 0.23 ps _1 eV otherwise, is plotted by open circles in Fig. 3. 
Even better agreement with the data is obtained from the TTM, assuming that the e-ph 
interaction strength F (e, e') entering Eq.Q smoothly vanishes as e and e' — > Ef, accounting 
for vf variation (and thus vf < v s condition) across the Fermi surface. This is implemented 
into the TTM simulation by approximating the factor D e (e) D e (e') F (e, e') in Eq. © with 
the symmetrized function (D e (e) D ie (e') + D e (e') D ie (e) — D ie (e) D ie (e')), where D ie (e) = 
D - £> exp[-(e/A') 2 ]. The resulting r ep (T), using A' = 24 meV, and goo = 4.5 x 10 14 
W/mol K is plotted by the solid line in Fig. 3. Indeed, extremely good agreement with the 
data is obtained, even though r r spans more than two orders of magnitude j^jj]. 

With the hypothesis that e-ph scattering is suppressed in the DOS peak, the experimen- 
tal observation of anomalous T-dependence of r T can be understood. Namely, at T < Tk 
C e (T) increases dramatically compared to normal metals. On the other hand e-ph relax- 
ation becomes more and more difficult as temperature is lowered since most of the electron 
relaxation should occur within the DOS peak, where the e-ph scattering is blocked by en- 
ergy and momentum conservation. Therefore, the thermalization between electrons and the 
lattice occurs very slowly, giving rise to the divergent r ep below Tk- While the presented 
model explains the main features of the data, i.e. the low-T divergence of r ep , there are still 
several issues requiring further experimental and theoretical effort. 

For example, in the simulations we considered a temperature independent peak in the 
DOS, assuming that many-body and correlation effects can be described by effective, T- 
independent model parameters. This may be an oversimplification of the physics of heavy- 
fermion systems. However, the relaxation time simulations and specific heat calculations of 
our phenomenological model depend only weakly on a T-dependent DOS, as long as the peak 
width in the DOS does not vary faster than temperature. Further, it would be interesting 
to investigate e-ph thermalization in Kondo insulators. Namely, due to the presence of the 
gap near E^ one would expect effects similar to the Rothwarf- Taylor bottleneck observed in 
superconductors Secondly, even more interesting effects are expected due to the strong 
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reduction of the screening at low frequencies (below the gap) which could lead to non- 
adiabatic phonons. 

In conclusion, we have utilized ultrafast optical spectroscopy to study the dynamics of 
photoexcited quasiparticles in HF compounds. We have observed a divergence in the e-ph 
thermalization time at low temperatures. We argue that the dramatic hundred-fold increase 
in the relaxation time at low temperatures in YbXCu 4 (and the lack of this quasi- divergence 
in the non-HF LuXCu 4 analogs) results from the largely increased DOS at Ef coupled with 
strongly suppressed scattering of heavy-electrons by phonons. 

We thank Kaden Hazzard for the specific heat data. This work was supported by US 
DOE. 

Figure Captions 
Figure 1 

Normalized photoinduced reflectivity data (solid symbols) on (a) LuAgCu 4 and (b) 
YbAgCu 4 at various temperatures, together with best fits (see text) to the data - solid 
lines. The data have been vertically shifted for clarity. 

Figure 2 

a) Temperature dependence of relaxation time r r on LuAgCu 4 (solid circles), together 
with the TTM prediction (dashed line) and the result of the numerical simulation (open 
circles) including the non-thermal electron distribution. Inset: C e (dashed) and C t (solid 
line) of LuAgCu 4 . 

Figure 3 

T-dependence of r r on YbAgCu 4 (solid circles), together with the corresponding r ep 
(multiplied by 10 for the presentation purpose) from simple TTM prediction (dashed line). 
Assuming suppressed scattering of heavy electrons by phonons due to vf < v s condition, 
very good agreement with the data is obtained: open circles present the result of numerical 
simulation, while solid line presents the TTM simulation - see text. Inset: experimentally 
determined C e (T) for YbAgCu 4 (open circles), together with calculated C e (T) based on the 
model D e (e) - solid line. C e (T) of LuAgCu 4 (dashed line) is added for comparison. 
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